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ABSTRACT 

We compare X-ray and weak-lensing masses for four galaxy clusters that comprise the top-ranked 
shear-selected cluster system in the Deep Lens Survey. The weak-lensing observations of this system, 
which is associated with A781, are from the Kitt Peak Mayall 4-m telescope, and the X-ray observations 
are from both Chandra and XMM-Newton. For a faithful comparison of masses, we adopt the same 
matter density profile for each method, which we choose to be an NFW profile. Since neither the X-ray 
nor weak-lensing data are deep enough to well constrain both the NFW scale radius and central density, 
we estimate the scale radius using a fitting function for the concentration derived from cosmological 
hydrodynamic simulations and an X-ray estimate of the mass assuming isothermality. We keep this 
scale radius in common for both X-ray and weak-lensing profiles, and fit for the central density, 
which scales linearly with mass. We find that for three of these clusters, there is agreement between 
X-ray and weak-lensing NFW central densities, and thus masses. For the other cluster, the X-ray 
central density is higher than that from weak-lensing by 2a. X-ray images suggest that this cluster 
may be undergoing a merger with a smaller cluster. This work serves as an additional step towards 
understanding the possible biases in X-ray and weak-lensing cluster mass estimation methods. Such 
understanding is vital to efforts to constrain cosmology using X-ray or weak-lensing cluster surveys 
to trace the growth of structure over cosmic time. 

Subject headings: cosmology: observations — galaxies: clusters: general — galaxies: clusters: indi- 
vidual (A781) — X-rays: galaxies: clusters — gravitational lensing 



1. INTRODUCTION 

Galaxy clusters have the potential to open a new win- 
dow on cosmology by serving as precision tracers of the 
growth of structure over cosmic time. The growth of 
structure can provide independent constraints on the 
matter density (fi M ), the dark energy density (Oa), and 
the dark energy equation of state (w), that would both 
verify our standard cosmological picture and take us fur- 
ther into unders t anding the nature of dark energy (e.g., 
iCarlstrom et al.l ((2002)). Utilizing galaxy clusters as 
tracers of structure growth largely relies on knowledge 
of cluster masses. Ideally cluster samples would have 
selection criteria based on mass, and mass estimates of 
clusters would be based on probes of their gravitational 
potential. However, most large samples of clusters that 
exist to date are selected on the basis of their trace 
baryons (i.e., visible light from galaxies or X-ray emis- 
sion from hot intracluster gas). Moreover, traditional 
probes of cluster mass (X-ray and optical) depend on 
the cluster's star formation history, baryon content, and 
assumptions about its dynamical state. Only recently 
have we obtained samples of clusters of significant size 
unbiased with respect to baryons and instead selected 
on the basis of their weak gravitational lensing shear. 

One such sample is provided by the Deep Lens Sur- 
vey (DLS), a deep BVRz' imaging survey of 20 square 



degrees (jWittman et al.1 l2002f ) . The observations were 
taken with the Cerro Tololo Blanco and Kitt Peak May- 
all 4-m telescopes. The primary goal of this survey is 
to study the growth of mass clustering over cosmic time 
using weak lensing. The DLS team has shown it is ca- 
pable of finding ne w galaxy clusters using the ir weak- 
lensing signal alone (jWittman et al.l 120011 . 120031 ). and it 
has presented its first sample of cluster candidates from 
the fi rst 8.6 square degrees of the survey (jWittman et al.1 
120061 ). The DLS survey should find ~40 clusters when 
completed. The CFHT Legacy Survey Deep has also pre- 
sented shear-selected clusters from a 4 square degree re- 
gion and the Garching-Bonn De ep Survey has presented a 
sample from 19 squa re degrees (jGavazzi fc Soucai]ll2007t 
ISchirmer etalll2007t l. 

We have been following-up a shear-ranked sample of 
DLS clusters with Chandra and XMM-Newton. One goal 
of this X-ray follow-up is to confirm that the DLS shear- 
selected cluster candidates are in fact true virialized col- 
lapsed structures. Preli minary analys i s for five of these 
clusters is presented in Hug hes et al.l (|2004T ). A further 
goal of this X-ray follow-up is to characterize the robust- 
ness of X-ray and weak-lensing cluster mass estimates 
and the biases inherent in X-ray and shear-selected sam- 
ples. This understanding is necessary in order to lay 
the groundwork for precision cosmology via larger X-ray 
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and weak-lensing cluster surveys (utilizing, for example, 
Constellation-X 1 and LSST 2 ). Such characterizations 
are facilitated by comparing weak-lensing mass estimates 
with X-ray mass estimates, as we elaborate on in §2. 

Below we report our weak-lensing and X-ray mass es- 
timates for our top ranked shear-selected cluster, A781, 
and three surrounding clusters. We discuss details of 
the weak-lensing and X-ray observations in §3, and the 
details of the weak-lensing and X-ray mass estimation 
methods in §4 and §5. In §6 we discuss our results, and 
in §7 we summarize our conclusions. 

2. BENEFITS OF INVESTIGATING THE RELATION 
BETWEEN X-RAY- AND WEAK-LENSING-DERIVED 
MASSES FOR SHEAR-SELECTED CLUSTERS 

An important issue for shear-selected clusters is pro- 
jection bias. The weak-lensing shear signal is sensitive to 
all the intervening matter between the background galax- 
ies and the observer. This leads to a possible projection 
bias of shear mass estimates as non-clust er line-of-sight 
matte r contaminates the sh e ar signal (e.g..lMetzler et all 
<fl999h : IWhite et all Ipoffi : Ide Putter fc White! (|2005fn 7 
X-ray observations provide an independent way to esti- 
mate the mass. Thus to quantify the extent of this con- 
tamination, we wish to compare X-ray mass estimates to 
weak-lensing mass estimates for clusters that are dynam- 
ically relaxed. X-ray observations are uniquely suited to 
this because they offer clear indications of a cluster's dy- 
namical state via X-ray images and temperature mea- 
surements. We impose the condition of relaxation be- 
cause X-ray mass estimates are based on an assumption 
of hydrostatic equilibrium and are likely invalid for highly 
unrelaxed systems. The comparison of X-ray to weak- 
lensing mass estimates will indicate how significantly pro- 
jection bias affects the latter. 

To study the bias (or absence of bias) inherent in X- 
ray and shear-selected cluster surveys, we first note that 
optical selection depends on star formation history and 
X-ray/Sunyaev-Zel'dovich selection depends on the heat- 
ing of the intracluster medium. It has been proposed 
that up to 20% of shear-selected clusters have not yet 
heated their intracluster m edium enough to be visible 
by cu rrent X-ray satellites (Weinber g fc Kamionkowskl 
2002). This is because a significant fraction of cluster- 
mass overdensities are likely nonvirialized and still in 
the process of gravitational collapse. These nonvirialized 
overdensities should produce much weaker X-ray emis- 
sion than that from a fully virialized cluster of the same 
mass. Differentiating this population from false-positive 
shear signals due to unrelated line-of-sight projections 
that appear as single larger mass concentrations, will be 
a challenge. Such a 'dark lens' cluster candidate was re- 
portedly found by lErben et al.l <|2000f) via weak-lensing 
observations ce ntered on Abell 19 42. This detection was 
followed-up by iGrav et "ail (|2001h in the infrared with 
no obvious luminous counterpart detect ed. Several more 
appar e nt 'dark lenses' are reported in iKoopmans et al. 
< 2000r> . lUmetsu fc Futamasd <|2000f l. and iMiralles et al. 
(2002). If such 'dark lenses' exist, there should exist 
a continuum of clusters between those which just satisfy 
Mxray < ^-weaklens ana ~ those which simply show no de- 
tectable X-ray counterpart to their weak-lensing signal. 

1 http:/ /constellation. gsfc.nasa.gov/ 

2 http:/ /www. lsst.org/lsst_homc.shtml 



Characterizing and quantifying the clusters for which 
Mxray < M wea ki ens will allow greater understanding of 
which clusters are missed by traditional samples and the 
percentage of false-positive detections that are inherent 
in shear surveys. 

The ratio of Mxray /^-weaklens ma y a -^ so P rove to be a 
good diagnostic of the dynamical relaxation of a cluster. 
Recent findings based on 22 high X-ray luminosity, low- 
rcdshift (0.05 < z < 0.31) clusters, selected on the basis of 
their high X-ray emission and targeted for weak-lensing 
follow-up with the ESO VLT, suggest X-ray cluster mass 
estimates larger than weak-lensing mass estimates pos- 
itively cojrelaAe_with clusters being dynamically unre- 
laxed (|Cvpriano et alj 12004). Naively one would expect 
this theoretically because events (such as mergers) that 
disrupt a cluster's equilibrium introduce transient shock 
heating of its intracluster gas. Calculating X-ray cluster 
masses using an assumption of hydrostatic equilibrium 
and a higher temperature than the cluster would have if 
relaxed, results in an overestimate of the true mass. How- 
ever, recent work based on hydrodynamic cluster simu- 
lations suggests X-ray mass estimates are biased low for 
unrelaxed clusters because only a portion of the kinetic 
energy of the merging system is converted into thermal 
energy of the intracluster medium, for even an advanced 
merger, while t he mass of the mergin g system has already 
increased (e.g.. lKravtsov et al. 2006). Comparing Mxray 
to M wea j-2 eng for our shear-selected clusters would de- 
termine whether X-ray mass estimates are biased high 
or low for unrelaxed clusters and whether this ratio can 
be used as a universal diagnostic of cluster dynamical 
state. Such a universal diagnostic would prove useful in 
investigating cluster evolution. 

Finally, there have been several reported instances of 
clusters that have an X-ray signal but no apparent weak- 
lensin g counterpart (|Cvpriano et alj 1200 4: Dahle et alj 
2002). We have detected such a cluster while following- 
up our highest shear-ranked cluster with XMM-Newton. 
This cluster did not appear in the original shear maps 
made for the DLS survey but is readily apparent in 
XMM-Newton observations. The inverse of 'dark lenses', 
negative weak-lensing detections are not unexpected 
since weak lensing is a less sensitive method of cluster 
searching as many galaxies need to be detected behind 
a cluster. Also mergers could potentially boost the X- 
ray signal of clusters otherwise below both current X-ray 
and weak-lensing thresholds. It is important for under- 
standing the limitations of weak-lensing surveys to ex- 
plore what is occurring in cases such as these. 

3. OBSERVATIONS 

3.1. Weak-Lensing Observations 

The Deep Lens Survey consists of five fields, each 
2° x 2° and isolated from each other. The two northern 
fields were observed using the Kitt Peak Mayall 4-m tele- 
scope, and the three southern fields were obtained with 
the Cerro Tololo Blanco 4-m telescope. Observing be- 
gan in November 1999 at Kitt Peak and in March 2000 
at Cerro Tololo. The deep BV Rz' images were ta ken 
with 8k x 8k Mosaic imagers (jMuller et alj I1998D on 
each telescope, which provided 35' x 35' fields of view 
with 0.26" pixels and minimal gaps between the CCD 
devices. The observing strategy was to require better 
than 0.9" seeing in the R band, so that this band would 
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IAU Designation 


Nickname 


CXOU J092026+302938 


Main Cluster 


CXOU J092053+302800 


Middle Cluster 


CXOU J092110+302751 


East Cluster 


CXOU J092011+302954 


Subcluster 


XMMU J091935+303155 


West Cluster 



TABLE 1 

IAU Designations for the clusters in the A781 cluster 
complex. 

have good, largely uniform resolution. When the seeing 
was worse than this, B, V, and z' images were taken. The 
source galaxy shapes were measured in the R band, and 
B, V, and z' image s provided color inform ation and pho- 
tometric redshifts. IWittman et al.1 (|2002l ) gives details of 
the fie ld selection and survey design, and Wittma n et al.l 
(2006) gives details regarding the image processing and 
convergence maps. 

A list of cluster candidates was compiled, based on the 
first 8.6 deg 2 of processed DLS data, and the candidates 
were ranked by their shear peak values. Multiple peaks 
within a 16' box were considered a single target for pur- 
poses of Chandra follow-up. A781 emerged as the top- 
ranked cluster candidate, with both DLS and archived 
Chandra observations indicating t hat this cluster was re - 
ally a complex of several clusters (|Wittman et al.l l2006). 
X-ray and optical follow-up of the A781 cluster complex 
was pursued as part of a larger follow-up program that 
will encompass a significant sample of DLS cluster can- 
didates. 

3.2. X-ray Observations 

We were awarded 15ks of XMM-Newton time in cy- 
cle 2 to get a closer look at our top-ranked DLS clus- 
ter complex. This observation took place on 04 April 
2003 (Obsid# 0150620201). In addition, Chandra had 
observed A781 on 03 October 2000 with the ACIS-I de- 
tector for a nominal exposure time of 10 ks (Obsid # 
534). The XMM-Newton and Chandra observations re- 
vealed that the A781 cluster complex consists of a large 
main cluster connected to a subcluster with two smaller 
clusters to its east and one to its west. We shall call 
the largest cluster the 'Main' cluster. The subcluster to 
its southwest appears in the act of merging with it (see 
Figure [l|. Just to the east of the Main cluster is another 
cluster, which we will refer to as the 'Middle' cluster, and 
within the same pointing, further to the east, is another 
cluster, hereafter 'East' cluster. The XMM-Newton ob- 
servation also presented us with a surprise. To the west 
of the Main cluster there appears to be one more cluster, 
which we will call the 'West' cluster. This cluster did 
not appear in the original DLS convergence maps made 
for the survey, and it is also, unfortunately, out of the 
field of view of the Chandra archive observations. Table 
[T] lists the IAU designations of these clusters. 

3.3. Optical Spectroscopy 

iGeller et al.1 (|2005l ) conducted a magnitude-limited (to 
B = 20.5) spectroscopic survey in this field. They re- 
port mean redshifts of 0.302, 0.291, and 0.427 for the 
Main, Middle, and Eas t clusters resp ectively (labeled as 
clusters A, B, and C in lGeller et aLT (120051) 1. These red- 
shifts were obtained from 163, 123, and 33 cluster mem- 



bers respectively. No redshift errors are quoted, however 
given the number of cluster members and the redshift 
errors in the individual galaxies, the systemic redshifts 
of the systems should be accurate to dz of 0.0002. The 
rest frame line-of-sight velocity dispersions were found to 
be a A = 674tg km s" 1 (Main), a B = 741±g km s" 1 

(Mid dle), and a c = 733±^ 2 km s" 1 (East) ([Geller et all 
2005). According to these velocity dispersions, these 
cluster components appear to be similar, a result we ex- 
amine further using our X-ray and weak-lensing data. 

The redsh i ft of the West cluster is not reported in 
IGeller et al.l |2005i). We o btained spectros copy of this 
cluster with Keck/LRIS (jOke et all |1995[ ) in longslit 
mode on 16 January 2007. We obtained secure red- 
shifts for two member galaxies, at a mean redshift of 
0.428 ± 0.001, though clearly the quoted error is itself 
highly uncertain with only two members. We also ob- 
served the East cluster in longslit mode on the same 
night, finding a mean redshift of 0.426 ± 0.003 base d 
on two members, in agreement with IGeller et al.1 (2005). 
Thus the East and West clusters are the only two com- 
ponents at the same redshift, but with a transverse sep- 
aration of 21' or 7.0 Mpc. Throughout this work we 
assume a A CDM cosmology of h = 0.71, £1a = 0.73, and 
Q M = 0.27 (|Spergel et alJl2003f ). 



4. EXTRACTED X-RAY TEMPERATURES AND GAS 
DENSITY PROFILES 

To obtain X-ray mass estimates of these clusters, we 
follow the standard practice of treating the intracluster 
gas as a hydrostatic fluid. This assumption is reasonable 
for dynamically relaxed clusters since collision times for 
ions and electrons in the hot gas are very short compared 
to times scales over which the gas heats or cools or the 
cluster gravitational potential varies. Assuming spherical 
symmetry, 

krT '(r) / din p(r) dhiT(r)' 



M(r) 



dim 



din 7 



(1) 



G/lTOp 

where /im p is the gas mean molecular weight, T(r) and 
p(r) are the gas temperature and density profiles, and 
M(r) is the total mass within a radius r. We assume the 
cluster gas follows a /3-model, 



pM=Po s [i+(fr 



(2) 



proposed bv lCavaliere fc Fusco-Femianol (|1978l ). In this 
model, the core radius, r c , is where the density is half the 
central density for a typical /3 of 2/3. To determine /3 and 
r c , we note that the X-ray emissivity is proportional to 
the square of the cluster gas density times the cooling 
function, i.e., e(r) cx A(T(r))(n e (r)) 2 . Integrating the 
emissivity through the cluster line of sight gives the X- 
ray surface brightness S. Assuming a p- model for the 
gas density and noting that the cooling function is close 
to constant over the range of typical cluster temperatures 
yields 
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e(r)dl cx 
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where I = yV 2 — b 2 and b is the projected radius (e.g., 
Sarazin (1986)). We fit the radial X-ray surface bright- 
ness profile with equation[3]to obtain estimates for (3 and 
r c (see §4.2). 
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The statistical quality of our data precludes the deter- 
mination of a radial temperature profile. In lieu of this, 
we assume an NFW profile for the cluster matter density 
given by 

Pm{t) = -, , — i—-^ (4) 
(r/r s )(l + r/r s y 

where pp is the central den sity and r s is the scale radius 
ijNavarro et alJll9M fl997h . This gives 

M (r) = 4 W f (ln(l + -) + 1 - l) (5) 
V r s 1 + rj r s J 

for the mass within a radius r. Using equation [T] we 
solve for the cluster temperature profile and then for 
the emission-weighted, projected, average temperature 
within a given aperture, which our data allow us to mea- 
sure relatively precisely. We compare this predicted tem- 
perature to our measured temperature determined from 
X-ray spectroscopy to find the best-fit value for po and 
thus the cluster mass. We describe the details of the 
X-ray spectroscopy below. 

4.1. X-ray Temperatures 
4.1.1. Analysis of Chandra Data 

We downloaded the Chandra data from the archive 
and used CIAO software tools for the initial data reduc- 
tion steps. The observation was carried out in full-frame 
timed exposure mode using very faint telemetry mode. 
The peak X-ray emission of the Main cluster was po- 
sitioned near the center of chip 13, some 3.5' from the 
on-axis "sweet" spot of the high resolution mirror. The 
merging subcluster to the west was imaged on chip 13, 
while the two other clusters toward the east were imaged 
on chip II. There was no cluster X-ray emission visible on 
the remaining two chips of the imaging array. We note 
that chip S2 of the spectroscopic array was also active 
during this observation, but we did not utilize these data. 
Event pulse heights were corrected for time-dependent 
gain, and all grades, other than 0, 2, 3, 4, and 6, were 
rejected. Information contained in the very faint mode 
data was used to reject non-X-ray background events. 
The light curve of the entire imaging array (minus obvi- 
ous cluster and unresolved emission) was examined and 
no time intervals of high or excessive background were 
found. The resulting live-time corrected exposure time 
was roughly 9900 sec. Figure Q] shows the 500-2000 eV 
band image of the ACIS-I data after exposure correction 
(for which the vignetting function was calculated at a 
monochromatic X-ray energy of 1 keV) and smoothing 
with a Gaussian of a ss 10". 

Spectral extraction regions for the clusters were de- 
termined (see Tables [5] and [3]) to optimize the signal to 
noise ratio of the resulting spectra. Annular regions sur- 
rounding each cluster were used to generate background 
spectra. Obvious point sources were excluded from both 
source and background regions, and cluster emission was 
excluded from all background regions. Weighted spec- 
tral response functions were generated for each source 
and matching background region, including instrumental 
absorption due to contamination build-up on the ACIS 
filters. 

Xspec (version 11.3) was used for the spectral anal- 
ysis. Fits were first done to the background spectra 



using a phenomenological model consisting of a non-X- 
ray background component (three gaussian lines and a 
power law to account for instrumental fluorescence lines 
and charged particles) and an astrophysical component 
(an absorbed power-law model to account for the unre- 
solved X-ray background). Inclusion of a soft thermal 
component (from nearby diffuse Galactic emission, for 
example) did not significantly improve the background 
fits, so it was not included. The absorption column den- 
sity was fixed to a value of Njj = 1.94 x 10 20 atoms 
cm~ 2 b ased on Galactic HI me asurements in this di- 
rection ijDickev fc Lockmanlll990h and the photon index 
of the astrophysical background component was fixed to 
r = 1.4. Only the normalization of this power-law model 
was allowed to vary. For the non- X-ray background, the 
gaussian line centroids and normalizations as well as the 
power-law index and normalization were allowed to vary 
freely. There were a total of nine free parameters for 
the fits to the background spectra. The source spectra 
included a redshifted thermal plasma model (mekal, in 
xspec parlance) to account for the cluster emission as 
well as the full component of background models just 
described. In all cases best fits were determined using 
the "c-stat" fit statistic, which is a likelihood figure of 
merit function appropriate for Poisson-distributcd data. 

Each pair of matched source and background spectra 
was fitted jointly with the background spectral compo- 
nents scaled between the source and background based 
on the ratio of exposure integrated over the extraction 
regions. In the joint fits the normalizations of the two 
background power laws plus the non- X-ray background 
photon index were allowed to vary. The other back- 
ground model parameter values were held fixed at values 
determined from the background fits alone. 

The best-fit temperature values and 1-a statistical un- 
certainties are given in Table El For the East and Middle 
clusters, the metal abundance was held fixed to 0.3 times 
solar. For the Main cluster, because of its higher statisti- 
cal level, this was allowed to vary yielding a best fit value 
of 0.27±0.15 times solar (1-er error). Redshifts were fixed 
to the values mentioned above. Figure [2] shows the best- 
fit spectra for these three clusters, where the dashed line 
represents the contribution from the background. 

4.1.2. Analysis of XMM-Newton Data 

The initial data reduction steps for the XMM-Newton 
data were completed using the SAS software tools. The 
data consists of observations from the three EPIC in- 
struments, MOS1, MOS2, and PN. To model the non-X- 
ray background from charged particles and instrumental 
fluorescence lines, we also obtained closed observations 
(observations with the filter wheel in the closed position) 
for each of the three EPIC instruments. Thus we were 
able to use an independent measure of the particle back- 
ground, differing from our Chandra analysis. 

Both our cluster data and the closed data were filtered 
by pattern and energy range. We kept single and double 
pixel events and events with pulse heights in the range 
of 300 to 12000 eV for MOS observations. For PN ob- 
servations, we only kept single pixel events and events 
with pulse heights between 300 and 15000 eV. We also 
chose the most conservative screening criteria (exclud- 
ing events next to edges of CCDs and next to bad pixels, 
etc.). Our cluster and closed EPIC data were also filtered 
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Chandra 


source 


background annulus (same center as source) 


Main Cluster 


09:20:24.8 +30:30:20.4 2.4' 


3.3' - 4.3' (minus Subcluster) 


Middle Cluster 


09:20:52.5 +30:28:08.4 1.6' 


2.3' - 3.3' (minus East) 


East Cluster 


09:21:10.9 +30:28:04.2 1.5' 


2.0' - 3.0' (minus Middle) 


Subcluster 


09:20:09.4 +30:30:02.5 0.9' 





TABLE 2 

chandra cluster and background annulus extraction regions. point sources excluded from the cluster regions are given 

in Table [3] 



R.A. 


Decl. 


Radius 


09:20:32.590 


+30:29:10.49 


4" 


09:20:29.774 


+30:28:55.73 


2" 


09:20:59.852 


+30:27:29.23 


8" 


09:20:50.765 


+30:29:22.97 


5" 


09:20:45.587 


+30:28:39.22 


8" 


09:21:06.817 


+30:27:40.20 


11" 


09:21:05.684 


+30:29:02.13 


6" 



TABLE 3 

Point sources excluded from the Chandra cluster regions. 



XMM-Newton 


region 


Main Cluster 


09:20:24.439 +30:30:21.12 2.5' 


Middle Cluster 


09:20:52.433 +30:28:12.74 1.4' 


East Cluster 


09:21:09.912 +30:27:58.31 1.1' 


Subcluster 


09:20:10.046 +30:29:57.17 1' 


West 


09:19:34.752 +30:32:00.88 1' 


Background Annulus 


09:20:12.521 +30:29:10.37 10.7' - 13.7' (minus East) 



TABLE 4 

xmm-newton cluster and background annulus extraction regions. point sources excluded from the cluster regions are 

given in Table [5] 



R.A. 


Decl. 


Radius 


09:20:33.209 


+30:28:58.48 


17" 


09:20:24.600 


+30:33:19.12 


24" 


09:20:21.809 


+30:30:35.14 


16" 


09:20:34.805 


+30:29:47.02 


12" 


09:20:25.531 


+30:33:39.11 


12" 


09:20:45.588 


+30:28:39.22 


12" 


09:21:06.816 


+30:27:40.20 


11" 



TABLE 5 

Point sources excluded from the XMM-Newton cluster regions. 



for soft solar proton flares. 

Figure Q] shows the 500-4500 ev band image of the 
EPIC data after exposure correction and smoothing with 
a Gaussian of a = 8". Extraction regions for the clus- 
ters and point sources were identified along with a back- 
ground annulus surrounding and away from any of the 
cluster regions (see Tables |3] and [S]) . Spectra were cre- 
ated of the cluster regions and background annulus for 
both the closed and cluster data, as were spectral re- 
sponse and effective area files. Resolved point sources 
were excluded from both the source and background re- 
gions, and cluster emission was also excluded from the 
background region. 

Xspec was used to fit the spectra of the background 
annuli in the closed and cluster data using a similar phe- 
nomenological model as described for the Chandra anal- 



ysis. The closed data background annulus was fit with 
several gaussians and three power laws to model the non- 
X-ray background. This best-fit model was used as a 
starting model to fit the spectrum of the cluster data 
background annulus, adding an absorbed power- law to 
model the unresolved X-ray background and an absorbed 
soft thermal (mekal) component to model the soft diffuse 
Galactic emission. The absorption column density and 
photon index of the astrophysical background model were 
fixed as mentioned above for the Chandra analysis, and 
the thermal component was given a fixed plasma temper- 
ature of 0.2 keV and a solar metal abundance. We linked 
the power-law norms of the particle background by the 
ratio of the power-law norms in the closed background 
data. This kept the slope and shape of the continuum 
fixed, but allowed the overall normalization to vary. The 



6 



Abell 781 XMM EPIC (0.5-4.5 klV) 




35' 



30' 



25' 



30° 20' 



Abell 781 Ctendra ACIS-KO.5-2.0 keV> 

m 




09 n 21 m s 



20 m 30 s 



00 s 



19 m 30 s 



R.A. (J2000) 



09 n 21 m s 



20 m 30 B 



R.A. (J2000) 



00 s 



Fig. 1. — Left panel: XMM-Newton 15 ks image of the A781 cluster complex. We refer to these clusters from left to right as East, 
Middle, Main, and West. The contours represent 1 X 10 -3 (white), 4.5 X 10 -4 , 2 X 10~ 4 , 8.9 X 10~ 5 , and 4 X 10 -5 cts/s/(4" square pixel). 
Right panel: Chandra 10 ks image of the A781 cluster complex. The contours represent 2.7 X 10 — 5 (white), 9.7 X 10 — 6 , and 3.5 X 10 — 6 
cts/s/(2" square pixel). Note the smaller field of view of the Chandra image, which only covers the three clusters to the east. 
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Fig. 2. — From top to bottom, spectra for the Main, Middle, 
and East clusters from Chandra. Solid lines represent the best- 
fit model, and dashed lines represent the contribution from the 
background. 

spectrum of the cluster data background annulus was 
fit by allowing the normalizations of the astrophysical 
background power-law and thermal component to vary 
as well as the normalization of the non- X-ray particle 
background. 

A best-fit joint model was created to fit the background 
annulus spectra for the three instruments simultaneously. 
The parameters of the astrophysical background model 
were kept in common between the instruments, but the 
parameters of the particle background differed. A sec- 
ond thermal component was added to the astrophysical 
background model to better fit the spectra at energies be- 



low the aluminum fluorescence line. The normalizations, 
temperatures, and abundances of the two thermal com- 
ponents were allowed to vary as well as the unresolved 
X-ray background power law. 

The source regions were fit drawing on the closed obser- 
vations to fit the non- X-ray background and the annulus 
spectra to fit the astrophysical background. The spectra 
of the cluster regions in the closed observations were fit 
using the closed background annulus best-fit model as 
a starting point. The spectra of the source regions in 
our data were fit using models starting where the parti- 
cle background normalization and the normalization of 
the strongest fluorescence line of each cluster region were 
set equal to those from the corresponding closed clus- 
ter region, scaled by the ratio of the normalizations be- 
tween the observed and closed background annulus. The 
normalizations of the weaker fluorescence lines were set 
equal to the normalization of the strongest line scaled by 
the ratio of the weak-line norm to the strong-line norm in 
the corresponding closed cluster region. The astrophys- 
ical background model was taken from the best joint-fit 
model for the background annulus, where the normaliza- 
tions of the background spectral components were scaled 
by the ratio of the exposure integrated over the source 
and background regions. The source spectra were mod- 
eled with a rcdshifted mekal thermal plasma model, with 
the abundance equal to 0.3 times solar and the temper- 
ature and normalization allowed to be free parameters. 

We fit the spectra of the source regions using these 
starting models by first fitting the normalization of the 
particle background using only events greater than 10 
keV, with the cluster model zeroed out. With this nor- 
malization frozen, we fit the temperature and norm of the 
cluster model using only events below 10 keV. The par- 
ticle background norm was then refit using events larger 
than 10 keV and this frozen cluster model. 
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Fig. 3. — Spectra for the Main cluster (top) and the East cluster 
(bottom) from the three XMM-Newton instruments. The Main 
cluster has the highest, and the East cluster the lowest, signal- 
to-noise of the four clusters in the pointing. Solid lines represent 
the best-fit model, and black, red, and green colors correspond to 
MOS1, MOS2, and PN instruments respectively. 
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Fig. 4. — Spectra for the Middle cluster (top) and the West 
cluster (bottom) from the three XMM-Newton instruments. Solid 
lines represent the best-fit model, and black, red, and green colors 
correspond to MOS1, MOS2, and PN instruments respectively. 



Having fit for the particle background of each source 
region for each instrument, the best-fit source models 
for the three instruments were combined to allow for a 
joint fit. Again for the joint fit, the astrophysical back- 
ground model was kept the same for each instrument 
except that the normalizations differed due to different 
exposure scalings. This fit was done using the best-fit 
particle background normalizations described above and 
those differing by ± 1-a for each instrument. In this 
way, we were able to model the systematic uncertainties 
arising from the particle background subtraction. 

The 1-a errors on the cluster temperature were ob- 
tained using a delta fit statistic of 1.0 for the one in- 
teresting parameter. The resulting best-fit temperatures 
and error bars are given in Table [6] We find good agree- 
ment between the Chandra and XMM-Newton best-fit 
temperatures given the statistical and systematic uncer- 
tainties. The XMM-Newton spectra for the four main 
clusters and the best joint-fit models are shown in Fig- 
ures [3] and [4] 



4.2. X-ray Surface Brightness Profiles 

A surface brightness profile was created by first deter- 
mining the surface brightness peak of each cluster. This 
was done with exposure corrected images in the 500-2000 
eV band for each instrument smoothed with a Gaussian 
of er = 8". The vignetting function for the exposure 
maps was calculated at the monochromatic X-ray energy 
of 1.25 keV for XMM-Newton and 1 keV for Chandra. 

The unsmoothed images and exposure maps from the 
three EPIC instruments and Chandra were used to gen- 
erate surface brightness profiles. Images and exposure 
maps from the separate EPIC focal plane detectors were 
summed to create a joint XMM image and correspond- 
ing exposure map. Exposure maps in all cases were made 
including the effective area, for proper weighting of the 
XMM detectors and for ease of comparing the XMM and 
Chandra profiles. 

We chose 40 radial bins of 8" each, centered around 
the surface brightness peak of each cluster to find profiles 
that extend out to 5.3'. In each radial bin, we summed 
the counts from the Chandra / 'joint XMM image and di- 
vided this sum by the total exposure in that bin from 
the corresponding exposure map. We calculated the er- 
ror on the surface brightness in each radial bin by using 
the small coun t statistic (1 + ^counts + 0.75)/exposure 
(|Gehreldll986l) . where units of exposure are given m sec 
cm 2 arcmin 2 . 

Using the radial bins farthest from the surface bright- 
ness peaks, we inferred the surface brightness due to the 
background. We then fit the surface brightness profiles 
to equation [3] by fixing the background and allowing f3, 
r c , and the overall normalization to vary as free param- 
eters. The profiles and best-fit models are displayed in 
Figures [5] and [6l Table [7] gives the (3 and r c best-fit val- 
ues for each cluster along with their 1-a statistical error 
bars. There is good agreement between the Chandra and 
XMM-Newton best-fit (3 and r c values given the statisti- 
cal uncertainties. 

5. MASS ESTIMATES 
5.1. X-ray Masses 

To determine each cluster's X-ray derived mass, we use 
equations [TJ [2j and [5] and our best-fit (3 and r c values to 
predict the cluster temperature profile for given values 
of the central density and scale radius in the NFW pro- 
file. The solution of the temperature profile requires a 
boundary condition. We assume the cluster temperature 
to be 1.25 keV at 3.5 Mpc for each cluster but allow this 
outer temperature and radius to vary generously over the 
range of 0.5 to 2 keV and 2 to 5 Mpc and fold this uncer- 
tainty into our error bars. We find that this variation in 
our boundary conditions contributes a negligible amount 
(less than 4%) to our error estimates. 

To compare each cluster's X-ray and weak-lensing- 
derived masses in an mutually consistent manner, we 
must assume the same matter density profile for each 
method. Since neither the X-ray nor the weak-lensing 
data is deep enough to well constrain both pg and r s 
in the NFW profile, we choose to estimate each cluster's 
concentration parameter, c cx l/r s , using results from hy- 
drodynamic simulations and an X-ray estimate of each 
cluster's mass assuming isothcrmality This mass esti- 
mate is accurate enough to give a reasonable estimate of 
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Fig. 5. — XMM-Newton surface brightness profiles for the Main (top, left), Middle (top, right), East (bottom, left), and West (bottom, 
right) clusters and best-fit models. The energy band used is 0.5-2 keV. 
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TABLE 6 

Integrated temperature estimates for the four clusters from fits to the XMM-Newton and Chandra spectra. Note the 

FIRST XMM-Newton ERROR GIVEN IS STATISTICAL AND THE SECOND IS SYSTEMATIC DUE TO BACKGROUND SUBTRACTION. 



Cluster 


XMM f} 


XMM r c (arcmin) 


Chandra P 


Chandra r c (arcmin) 


East Cluster 


0.811H 


i - ly -0.24 


0.681S;™ 


0.94^ 


Middle Cluster 


51+ lulb 


U -' 1 -0.1B 




0.99^ 


Main Cluster 


0.871^ 


i.82+;; ; - 


0.88t u; ^ 




West Cluster 




u - J - L -0.03 







TABLE 7 

Best-fit P and r c values for the four clusters from fitting XMM-Newton and Chandra surface brightness profiles assuming 

a /3-model for the gas density. 
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Fig. 6. — Chandra surface brightness profiles for the Main (top, left), Middle (top, right), and East (bottom, left) clusters and best-fit 
models. The energy band used is 0.5-2 keV. 
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1.48 
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TABLE 8 

Best-fit Chandra and XMM-Newton mass estimates within r 5 oo assuming an isothermal /3-model. 
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the concentration parameter since the concentration is 
a slowly varying function of cluster mass. The concen- 
tration parameter is here defined as r5oo/r s , where rsoo 
is the radius within which the cluster mean density is 
500 times the critical density. We again use the isother- 
mal case to estimate rsoo, which yields an estimate of 
r s . Since we are primarily interested in the comparison 
between X-ray-derived and weak-lensing-derived masses, 
the accuracy of the scale radius is less important than 
the fact that we used the same scale radius in deriving 
the masses using both methods. A change in the scale 
radius may introduce a systematic shift in the masses 
but will not alter significantly their ratio. 

5.1.1. Masses Assuming an Isothermal-[3 Model 

We determine X-ray isothermal /3-model mass esti- 
mates for these clusters as an intermediate step to es- 
timate the NFW scale radius for each cluster. If we as- 
sume that each cluster is isothermal with a gas density 
given by equation [2j then equation [1] can be rewritten to 
give 



M{r) = 3/3 



kTr 



Gpm p 
1.13 x 10 15 /3 



(r/r c 



1 + (r/r c ) 
T x r 



10 keV Mpc 



( r / r cY 



l+(r/r c y 



M Q , (6) 



wher e we set p = 0.59 for the cluster gas (|Evrard et al.1 
1996). Using the best-fit /3, r c , and Tx for each cluster 
from Chandra and XMM-Newton observations, we calcu- 
late M500 for each cluster, which is the mass within rsoo- 
The best-fit rsoo and M500 for each cluster from each X- 
ray satellite are listed in Table [H This table is without 
error bars since we only use these values to estimate a 
reasonable r s for each cluster. 

We use the above XMM r 500 and M 50 o values and 
predictions for the concentration as a function of mass 
and red shift derived from hydrodynamic cluster simu- 
lations (jDolag et al.ll2004T ) to obtain an est imate of r s . 
The c (M, z) relation used here is from the iDolag et al.l 
( 2004) ACDM cosmological simulation, and we converted 
our masses and radii to the definitions used in that work 
to determine r s . This concentration relation has a rea- 
sonable agreement with Chand ra observations of nearby 
clusters (|Vikhlinin et al.ll2006l ). We choose to use XMM 
/3-model masses since we have XMM observations for all 
four clusters. Table [5] lists the best estimates of r s . We 
hold this scale radius fixed and keep it in common when 
deriving masses from weak-lensing and X-ray methods 
and focus on the relative difference in mass estimates 
using these two methods. 

5.1.2. Masses Assuming an NFW Profile 

Given an estimate of r s , /3, and r c for each cluster, we 
vary po and calculate the emission-weighted, projected, 
average temperature within an aperture (as discussed in 
2]). The apertures for each cluster and instrument are 
given in Tables O and 2] We match the predicted temper- 
atures to our measured temperatures in Table [5] to find 
the best-fit central densities. The 1-a errors on po are 
calculated by varying the measured Tx, P, and r c param- 
eters within their 1-a error bars and treating /3 and r c as 



correlated variables. The variation in the boundary con- 
dition introduces a negligible error on po as mentioned 
above. Best-fit po measurements and 1-a errors for Chan- 
dra and XMM observations are given in Table [9] There 
is good agreement between the two instruments, and we 
calculate a weighted average to give a combined X-ray po 
for each cluster, which we also list in Tabled] Given p 
and r s , X-ray masses are calculated using equation [5] and 
listed in Table [TU1 These masses are calculated using the 
combined X-ray po for all the clusters except the West 
cluster, for which only the XMM po is available. 

5.2. Weak-Lensing Masses 

The weak lensing mass estimates are based on the full 
DLS exposure time of 18 ks in R and 12 ks in BVz' , 
rather than the partially complete im aging used for clus- 
ter discovery in lWittman et all |2006). Otherwise, image 
processing was as described in detail in that paper. The 
R data were taken in good seeing conditions (FWHM 
< 0.9") and are used to measure the galaxy shapes. The 
A781 complex spans two contiguous DLS "subfields" or 
pointing centers. The FWHM of the final R images, af- 
ter circularizing the PSF and co-adding 20 exposures, 
is 0.78" and 0.74" for the two subfields involved. The 
BVz 1 data are used only to provide color information for 
photometric redshifts. 

We ext racted shear catalog s using a partial implemen- 
tation of iBernstein &: Jarvisl (|2002f ). T his implementa- 
tion a ppeared as the "VM" method in [Hcvman s et al.l 



(2006), which compared different weak-lensing methods 
on a set of simulated sheared images. After correcting for 
stellar contamination which was present in that dataset 
but not here, the VM method yielded 89% of the true 
shear in those images. In this work, we divide the VM re- 
sults by 0.89 to more closely approximate the true shear. 
However, we recognize that this correction factor is likely 
to be data-dependent, and we therefore assign a system- 
atic error of 10% to the shear calibration. 

We derived photometric redshifts using BPZ (|Bemtezl 
[2001 with the HDF prior. We optimized the t emplates 
using a subset of the SHELS dGeller et al.ll2005h spectro - 
scopic sample and the procedure of lllbert et al.1 (|2006h . 
Complete details are discussed elsewhere (Margoniner et 
al., in preparation). To assess the accuracy of these pho- 
tometric redshifts, we turn to an independent spectro- 
scopic sample, consisting of all redshifts in DLS field F2 
in the literature, as tabulated by the NASA/IPAC Ex- 
tragalactic Database. We find 328 galaxies with spec- 
troscopic redshifts that match our cleaned photometry 
(unsaturated, not near bright stars, etc), spanning the 
redshift range 0.02 — 0.70. The resulting rms photomet- 
ric redshift error per galaxy is 0.047(1 + z), with a bias 
of —0.017(1 + z), and no catastrophic outliers. 

Because shear is nonlocal and mass maps tend to be 
highly smoothed, the presence of one clump may affect 
the appar ent mass density of a nother clump. This is 
true of the lWittman et all (|2006f ) maps, and we eliminate 
that effect here by simultaneously fitting axisymmetric 
NFW profiles to the four X-ray positions. The model 
fitting takes into account the full three-dimensional posi- 
tion (RA, DEC, z) of each source galaxy. The per-galaxy 
imprecision in z is not important because the lensing 
kernel is very broad, and because each galaxy is a very 
noisy estimator of the shear: a 0.1 shift in source red- 
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shift changes the modeled shear by much less than the 
per-galaxy shear error. For each NFW model, RA and 
DEC were fixed by the X-ray position, z was fixed by the 
spectroscopy, and r s was fixed to the value used for the 
X-ray fitting . Thus only one pa r amete r, po, was fit for 
each model. IWright &: Brainerdl <|2000l ) give expressions 
for shear induced by an NFW profile. 

Of the ~350,000 galaxies in the 2°x2° DLS field con- 
taining A781, we limited the fit to galaxies within 15' 
of a clump center, for computational efficiency and be- 
cause more distant galaxies may be influenced more by 
other clusters than by those in the A781 complex. We 
also cut on photometric redshift, because cluster mem- 
bers scattering to higher redshift would reduce the esti- 
mated shear (slightly, because of their low inferred dis- 
tance ratio), while cluster members scattering to lower 
redshift would have no effect. We did fits with cuts 
at Zphot > 0.35 (just behind the richer, lower-redshift 
clumps) and z p hot > 0.62 (3er beyond the higher-redshift 
clumps), and found a difference of < 0.2a in the fit- 
ted parameters. The lenient (strict) cut yielded 30137 
(22173) galaxies, or 23 (17) arcmin -2 over 1320 arcmin 2 , 
although ~10% of this area was masked due to bright 
stars. We adopt the strict cut to avoid any question of 
contamination. The resulting source catalogs show no 
increased density near the clump centers. 

Because shear from an NFW profile is linear in po, we 
used singular valu e decomposition (SVD) as described in 
iPress et al.1 ((T992)) . This solves the general linear least 
squares problem in one pass, with no iteration required. 
The galaxies were given equal weights in the fit, because 
the VM shear method does not assign weights to galaxies. 
However, the importance of a galaxy in determining the 
fit still depends on its position and redshift, through the 
model's dependence on position and redshift. We then 
corrected for the fact that the observable in weak lens- 
ing is not the shear 7, but the reduced shear 7/(1 — k) 
(where k is the convergence) as follows. We computed 
the convergence of the best-fit model at the location of 
each source galaxy, constructed a reduced-shear model, 
redid the linear fit, and iterated. In the first iteration, 
this resulted in a ~5% correction to the fit parameters. 
In the second iteration, the correction was only ~0.2%, 
much smaller than the fitted parameter statistical uncer- 
tainties, and therefore the reduced-shear fit was deemed 
to have converged. 

The fitted po 's and their uncertainties are listed in Ta- 
ble [9] The uncertainties output by the SVD routine were 
confirmed by 1000 bootstrap resampling realizations. Al- 
though the fit is not x 2 -driven, we can define a A% 
statistic to see how much of the variance in the data is 
explained by the model. There are 44326 degrees of free- 
dom because each galaxy has two shear components. For 
each component of shear, we compute the rms without 
any fit and define that as the uncertainty associated with 
each galaxy. This results in an initial \ 2 °f 44326, which 
is decreased by 40.4 when the four-parameter model is 
subtracted off. The chance of this happening randomly is 
< 10~ 7 for Gaussian distributions. The remaining vari- 
ance is due mostly to the intrinsically random distribu- 
tion of galaxy shapes (shape noise) and shape measure- 
ment errors, although a small amount may be attributed 
to additional structure in the field as indicated by the 
following mass maps. 




09 h 21 m ' " 20 m 

Fig. 7. — Lensing data, model, and residuals. Top: mass map 
of the area. Grayscale and contours are in arbitrary units, but 
the same units are used for all the panels. Middle: best-fit lensing 
mass model. Bottom: mass map made from residual shear after 
subtracting off the quadruple NFW profile fit shown in the middle 
panel. 

We show the fit visually in Figure [7] The top panel 
shows a mass map ma de using the me thod described in 
IWittman et all (|2006h (originally from iFischer fc Tvsonl 
( 1997I )). the middle panel shows a similar map made 
from the model shears, and the bottom panel shows a 
map made from the residual shear after subtracting off 
the fit. The Main, Middle, and East clumps have been 
mostly subtracted, but the West clump has not been well 
modeled. There is also an unmodclcd mass clump just 
northwest of the East clump (RA w 9:21, Dec w 30:33), 
which appears to have some associated galaxies. We cau- 
tion that the mass map is a sanity check rather than a 
quantitative indicator of goodness-of-fit, because it does 
not fold in source redshift information. 

In Table O weak-lensing errors are of two types, the 
first is statistical and the second is systematic. For the 
statistical errors, we performed bootstrap resampling to 
estimate the covariance of the cluster mass estimates. 
The masses of neighboring clusters (in projection) are 
anticorrelated, because the observed shear in a region 
must include the sum of the model shears from the neigh- 
bors. The errors given here for each cluster are after 
marginalizing over the allowed values for the other clus- 
ters. Therefore, the error on the total mass of all the 
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Cluster 


Estimated t s 
(Mpc) 


XMM po 
(10- 26 g/cm 3 ) 
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3.2 ± 0.2 
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TABLE 9 

Estimated scale radii of NFW mass profiles and best-fit NFW central densities from Chandra, XMM- Newton, and DLS. 
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X-ray r 50 o 
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Weak-lensing rsoo 
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Weak-lensing M500 
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1 - ii -0.5 
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TABLE 10 

X-RAY AND SHEAR MASSES WITHIN rsoo FOR THE FOUR CLUSTERS ASSUMING AN NFW MATTER DENSITY PROFILE. 



clusters would be smaller than the quadrature sum of 
the errors given here. Also, the covariance will affect the 
comparison of the lensing p$ for one of the four clusters 
to that for any of the others. 

Systematic errors include shear calibration, source 
rcdshift calibration, mass sheet degeneracy, and resid- 
ual cluster member contamination. We assign a shear 
calibration systematic uncertainty of 10% as explained 
above. We explored the effects of source redshift errors 
by fitting to different source redshift ranges. We find 
that only if we include photometric redshifts beyond 1.4, 
where the 4000 A break redshifts out of the z 1 filter, do 
the results change by as much as 10%. Even in that case, 
the results appear to change randomly rather than sys- 
tematically, with some clumps increasing in fitted mass 
and others decreasing. Given these results, and the small 
bias in photometric redshifts when compared to the spec- 
troscopic sample from the literature, a systematic error 
of 5% in mass due to redshift calibration errors is very 
conservative. 

We estimate residual cluster member contamination by 
examining the source galaxy areal density around the 
richest (Main) clump. There is an excess in the cen- 
tral 2' radius. Comparing the redshift distribution in 
that area to a control annulus, we find 58 excess galaxies 
at redshifts 0.35-0.42, presumably cluster members with 
> la- photometric redshift errors. In the fit, all galax- 
ies were weighted equally, so their relative importance is 
determined by their proximity (in three dimensions) to 
where the model reduced shear is large. These galaxies 
are near the projected center of the lens, but at low in- 
ferred distance ratio. We compute their effective weight 
as the square of the model reduced shear, and compare 
their total weight to the rest of the galaxies within 5' of 
the Main clump. The summed weight of the interloping 
galaxies is only 0.010 times that of the valid source galax- 
ies, which presumably resulted in a 1% underestimate of 
the mass, much smaller than the other systematics. 

The last systematic involves mass sheet degeneracy 
and projection of unrelated structures near the line of 
sight. // our assumption of an NFW profile is correct, 
then mass sheet degeneracy is not an issue, because we 



are fitting the profile rather than empirically determin- 
ing departures from a baseline. Measuring shear to as 
large a radius as possible would help check the profile 
assumption, just as it would help reduce mass sheet de- 
generacy. However, going to large radii increases the 
chance of including so me unrelated st r uctur e projected 
near the line of sight. iMetzler et ail (|200lD (hereafter 
M01) characterized this by measuring lensing masses of 
clusters in large-scale numerical simulations as a function 
of data set radius. The truncation radius used in the fit- 
ting here, 30', corresponds to a transverse separation of 
about 5 Mpc, for which M01 find a scatter of about 7%. 
(Note that MOl's systematic offset in the population is 
not relevant here because profile fitting is less biased than 
their aperture densitometry.) Scatter in the population 
becomes a systematic when considering a single cluster. 
However, by simultaneously fitting multiple clumps, this 
systematic is probably already greatly reduced. We em- 
pirically test this effect by varying the truncation radius 
and find that the results can change by up to 5%. We 
therefore assign a systematic of 5% due to this effect. 
The larger systematic is likely to be in the profile as- 
sumption. This systematic is likely to be on the same 
order as the mass sheet degeneracy systematic incurred 
if the pro file assumption were dropped, which is quoted 
as 20% by iBradac et all (|2004h . 

In summary, the systematics include 20% for the pro- 
file assumption and/or mass sheet degeneracy, 10% for 
shear calibration, 5% for projections, and 5% for photo- 
metric redshift errors. We include only the latter three 
systematic errors in Table |H1 where we assume an NFW 
profile accurately describes the matter density. This al- 
lows us to compare the X-ray and weak-lensing po values. 
In table [TOj we present a summary of the r^oo and M500 
values, where the the weak-lensing statistical and sys- 
tematic error bars on po are added in quadrature. It 
should be kept in mind that the absolute masses of the 
clusters depend on the profile assumption. 

6. DISCUSSION 

Comparison of the X-ray and weak-lensing po values 
indicates that these values are in agreement for the East, 
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Middle, and West clusters. This agreement suggests that 
line-of-sight mass contributions have not significantly bi- 
ased the weak-lensing measurements. 

For the Main cluster, the X-ray-derived po is higher 
than that from weak-lensing by about 2a. The X-ray 
images of the Main cluster suggest it may be undergoing 
a merger with a subcluster (see Figured]), and thus may 
be out of hydrostatic equilibrium. From the literature it 
is not obvious whether cluster mergers are expected to 
bias X-ray mass estimates high or low. Weak-lensing ob- 
servations of 22 X-ray bright clusters at 0.05 < z < 0.31 
found X-ray temperatures higher by 40 — 75% than those 
inferred from weak-len sing, for clusters with Tx > 8 keV 
(jCvpriano et al.ll2004l ). The largest discrepancy in this 
sample was for the two clusters with the highest X-ray 
temperatures (Tx ~ 13 keV), which both show signs of 
being out of dynamical equilibrium. It is a reasonable ex- 
trapolation to presume that all the clusters with Tx > 8 
keV in their sample are unrelaxed, with the temperature 
of their intracluster medium boosted by shocks due to 
in-falling groups and mergers with other clusters. 

However, recent hydrodynamic simulations indicate 
that unrelaxed cluster temperatures should be lower 
than those for relaxed clusters of th e same mass 
(jKravtsov et alj 120061: iNagai et alj 120071) . The reason- 
ing for this is that over the course of a merger, the 
mass of the system increases faster than the conversion 
of the kinetic energy of the merging sys tems into ther- 
mal e nergy of the intracluster medium (jKravtsov et alJ 
2006). It has also been suggested from hydrodynamic 
simulations of merging clusters that X-ray mass esti- 
mates based on hydrostatic equilibrium can be biased 
high close to core-crossing, where a temperature boost 
occurs due to shocks, but can be bi ased low just be- 
fore and after this temperature spike (|Poole et al.l 120061 
iPuchwein &; Bartelmannll2007| ). This latter scenario is 
possibly supported by recent X-ray and weak-lensing 
observations of 10 X-ray luminous clus ters at z ~ 0.2 
(|Zhang et alj|2007l iBardeau et al.ll2007l) . The X-ray ob- 
servations of this sample were carried out with XMM- 
Newton, and the weak-lensing is from ground-based 
imaging with the CFH12k camera on CFHT. These au- 
thors report that four out of these ten clusters have con- 
sistent X-ray and weak-lensing mass estimates. Of the six 
clusters that show a mass discrepancy, half have higher 
and half have lower X-ray as compared to weak-lensing 
masses. 

It is possible that the Main cluster in the A781 clus- 
ter complex is close to core-crossing as our X-ray mass 
estimate is biased high. It is also possi ble that a selec- 
tion e ffect is occurring in the sample of ICvpriano et aTl 
(|2004h where the clusters with T x > 8 keV are m a 
state of high temperature boosting and thus close to core- 
crossing. Deeper X-ray observations of this A781 cluster 
and further comparisons of weak-lensing and X-ray mass 
estimates of known merging clusters will help to clarify 
the biases expected from dynamically unrelaxed systems. 

The West cluster did not appear in our original shear 
maps, and we confirm here that it is a low significance 
weak-lensing detection (1-2<t). It was detected in the X- 
ray by chance as it fell within the same XMM-Newton 
pointing as the other three clusters. However, weak- 
lensing and X-ray data yield consistent mass estimates. 
We measure a very small core radius for this cluster 



(r c = 0.31'), consistent with most of the emission com- 
ing from a compact core. Northwest of the East cluster, 
we also find an enhancement in both the galaxy distri- 
bution and lensing signal, and we find some indication of 
X-ray emission from that region. Further X-ray observa- 
tions would allow us to study these lower mass systems 
in more detail. 

Based on the limited information available, the ve- 
locity dispersion values appear broadly consistent with 
the East , Middle, and Mai n clusters having nearly equal 
masses dGeller et all 120051) . This agrees with the simi- 
larity in weak-lensing masses we present here (see Table 
[TO]) . A detailed comparison of the cluster masses de- 
rived above with the velocity distributions is beyond the 
scope of this work. Study of the velocity distribution of 
the potentially merging component (Main cluster) could 
shed light on the ep och and geometry of the merger (see 
iGomez et al.1 (|200a) V 

7. CONCLUSION 

Many cluster surveys will take place over the next few 
years at several different wavebands. Already consid- 
erable samples of X-ray and optically selected clusters 
have been compiled. In addition, sizeable microwave and 
shear-selected samples are close at hand. These surveys 
can probe with precision the growth of structure over 
cosmic time, and thereby open a new window on cos- 
mology. However, the two main hurdles to overcome are 
relating cluster observables to mass and characterizing 
the sample selection. By comparing weak-lensing and X- 
ray observations and mass estimates for clusters in the 
DLS shear-selected sample we hope to understand the 
systematic biases in both mass estimation methods and 
modes of cluster selection. 

An analysis of the top shear-ranked mass distribution 
in the DLS sample reveals a complex of four clusters, 
the largest of which can be identified as A781. The 
four clusters are at distinctly different redshifts, as de- 
termined by optical spectroscopy, and the X-ray images 
suggest three are dynamically relaxed while the largest 
cluster appears to be merging with a small subcluster. 
Masses from both X-ray and weak-lensing observations 
were determined assuming an NFW profile for the mat- 
ter density. Since neither sets of observations were deep 
enough to well constrain both the central density and 
scale radius of each cluster NFW profile, we estimated 
the scale radii using X-ray-derived isothermal /3-model 
mass estimates and a relation describing concentration as 
a function of mass and redshift derived from cosmological 
hydrodynamic simulations. For each cluster profile, the 
same scale radius was used to determine the best-fit X- 
ray and weak-lensing central densities. We focus on the 
difference in central densities derived with each method 
as the central density scales linearly with mass. 

We find that three out of the four clusters show agree- 
ment between their X-ray and weak-lensing derived cen- 
tral densities. The fourth and largest cluster has an 
X-ray derived central density higher than that derived 
from weak-lensing by about 2a. This discrepancy is most 
likely due to the cluster's disrupted dynamical state. Re- 
cent weak-lensing observations of X-ray selected clus- 
ters and hydrodynamic simulations leave some ambiguity 
about whether dynamical disruption via cluster merg- 
ers biases X-ray mass estimates high or low. The di- 
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rection of the bias may be related to the stage of the 
merger, e.g., whether it is close to core-crossing or more 
advanced. Deeper X-ray observations of this cluster to 
better resolve the merger and further comparisons be- 
tween weak-lensing and X-ray-derived masses of known 
merging clusters will shed greater light on this issue. 

Initial steps are being made by many groups to over- 
come the above mentioned hurdles regarding the use of 
clusters as precision tracers of structure growth. Our col- 
laboration, for example, has Chandra and XMM-Newton 
data on a larger sample of DLS shear-selected clusters 
that we will be reporting on in future publications. Hope- 
fully all these efforts will open a new window through 
which we can understand our universe. 
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